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Abstract 

A rephasing invariant parametrization is introduced for three flavor neutrino mixing. For neu- 
trino propagation in matter, these parameters are shown to obey evolution equations as functions 
of the induced neutrino mass. These equations are found to preserve (approximately) some char- 
acteristic features of the mixing matrix, resulting in solutions which exhibit striking patterns as 
the induced mass varies. The approximate solutions are compared to numerical integrations and 
found to be quite accurate. 
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I. INTRODUCTION 



Flavor mixing plays a central role in the physics of flavors. For quarks, the CKM (Vckm) 
matrix has stood the test of time and is found to be sufficient in describing all of the relevant 
physics. Similarly, the PMNS (V u ) matrix has been used to analyze neutrino oscillation with 
no known discrepancies. Mathematically, both matrices belong to elements of U(3), the 
3x3 unitary matrices. Physically, for quarks, since the phases of individual quark states 
are unobservable, the rephasing transformation, Vckm PVckmP', where P and P' are 
arbitrary phase matrices, leaves the physics unchanged. Thus, only the rephasing invariant 
part of Vckm is physical. For V v , while the charged lepton phases are unobservable, for 
Majorana neutrinos actually there are two observable, CP- violating, phases However, 
as long as one restricts oneself to lepton number conserving processes, such as in neutrino 
oscillations, these phases also become unphysical so that the physical V v is again rephasing 
invariant. 

Related to the rephasing invariance of mixing matrices is their parametrization. While 
it may appear that the choice of parametrization is not important, since, at the end of the 
day, the physical quantities must be grouped into rephasing invariant combinations. How- 
ever, when one deals with a situation where relations between parameters are considered, a 
particular choice may be advantageous over others. For instance, when the mixing depends 
on the energy scale, as in the RGE for mass matrices, we have a set of evolution equations 
relating parameters at neighboring scales. Another example deals with neutrino mixing in 
matter. Here, the mixing depends on the density of matter and neutrino energy, contribut- 
ing to an induced neutrino mass. One can establish the relation between parameters for 
neighboring densities, resulting in a set of evolution equations as a function of the neutrino 
effective mass. They are very similar to those of the energy scale as described by a set of 
RGE. It turns out that, in both cases, the use of explicitly rephasing invariant parameters 
simplifies the evolution equations. 

In the following we will derive a set of evolution equations, as a function of the effective 
mass of neutrinos, for neutrino parameters in matter. These equations are based on the 
use of rephasing invariant parameters developed earlier. We find that they have simple, 
analytic, albeit approximate, solutions. It is interesting that the parameters in matter 
preserve a number of salient features of those in vacuum, resulting in a matter-dependent 
PMNS matrix that can be grasped at a glance. 

The paper is organized as follows. Section II is a brief summary of the rephasing invariant 
parametrization that is adopted in this work. In stead of directly solving the eigenvalue 
problem, we derive in Section III the evolution equations for the neutrino mixing parameters 
and masses from the effective Hamiltonian in matter. Certain well-known invariants are 
also derived using the symmetric properties of the equations. Section IV is devoted first 
to solving the two-flavor problem using this rephasing invariant formulation, and then the 
three-flavor case. Making use of the known properties of measured neutrino parameters, 
analytic, approximate, solutions are obtained. In Section V, the accuracy of the solutions 
are confirmed by comparison with numerical integration of the equations. Section VI is 
the summary. In appendix A, we also derive the neutrino transition probabilities in matter 
using the adopted rephasing invariant parametrization. 
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II. REPHASING INVARIANT PARAMETRIZATION 



In this section, we briefly summarize the rephasing invariant parametrization introduced 
earlier for quark mixing [2], and will now be adopted for neutrino mixing, valid for lepton 
number conserving processes. 

For the PMNS matrix (V), without loss of generality, we can impose the condition detV = 
+ 1. There are then a set of rephasing invariants 

r^fc = ViiV2jVsk = Rijk — iJ, (1) 

where their common imaginary part can be identified with the Jarlskog invariant J ^ . Their 
real parts are defined as 



(#123, #231, #312; #132, #213, #32l) = %2, X 3 ] Vl, 2/2, Hz) ■ 



(2) 



These variables are bounded by ±1: —1 < (xi,yj) < +1, with yj < X{ for any They 
satisfy two constraints 



detV = (xi +x 2 + x 3 ) - (yt + y 2 + 2/3) : 
{x x x 2 + x 2 x 3 + £3X1) - (2/12/2 + 2/22/3 + 2/32/1) 

In addition, it is found that 

J 2 = xix 2 x 3 - 2/12/22/3- 
The (x,y) parameters are related to \Vij\ 2 by 



0. 



(3) 
(4) 

(5) 
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(6) 



One can readily obtain the parameters (x, y) from W by computing its cofactors, which form 
the matrix w with w T W = (detW)I, and is given by 



w 



f x x + yi x 2 + y 2 x 3 + y 3 ^ 

x 3 + 2/2 xi + 2/3 x 2 + yi 
\x 2 + y 3 x 3 + yi xi+y 2 J 



(7) 



Eqs. ([6]) and ((Zj) establish the close relationship between the two rephasing invariant 
parametrizations (x,y) and \Vij\ 2 . Besides the obvious difference in the number of con- 
straints (two for (x,y) and five for |V^| 2 ), the set (x,y) has built-in symmetry amongst the 
three states considered, which, as we will see, helps to make the evolution equations simpler. 
For the PMNS matrix in vacuum, its elements squared are well-approximated by 



Wn 



\ 



2 



(8) 



with (e, 77, <C 1. Wq reduces to the tri-bimaximal Q matrix when e = rj = (3 = 
£ = 0. If we allow the parameters (e, 77, /?,£) to take on arbitrary values, the matrix above 
can be used as a general parametrization of the mixing matrix. Also, it is related to the 
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familiar "standard parametrization" jEf by Sf 3 = e 2 , S\ 2 — § + 1^2 > ^23 = I + ll^?; anc ^ 
2/3 = (S 2 23 - C 2 2 3 )[-|C 2 3 + C\ 2 - S 2 12 Sf 3 ] + 45 12 C 12 5 13 C 23 5 2 3Cos07so that, if (e,rj,i) « 1, 

The matrix PFo m Eq. ([8]) exhibit several interesting features. When e = ^ = /3 = ^ = 0, 
we find £10 = 1/3, x 2 o = 1/6, £30 = 0, and + yio = (z = 1,2,3). The conditions 
£30 = 2/30 = come from W13 = (so also V13 = 0). The conditions + Uio = are 
equivalent to W 2 i = W^i From known experimental bounds, for non- vanishing (e, rj,/3,£), 
these conditions are valid to C?(10~ 2 ). 



III. EVOLUTION OF NEUTRINO MIXING PARAMETERS 



It is well-established that neutrino mixing is modified by the presence of matter Q. Their 
effect has been used in the analyses of solar neutrinos, and is expected to impact those of 
the supernova neutrinos, when and if they become available. Closer to home, there is a 
plethora of long baseline experiments either in operation or in the planning stage. For these 
studies, it is essential to include the matter effects in order to understand neutrino mixing 
at the fundamental level. 

In the literature, effort has been devoted to solving problems along this line (see, e.g., 0). 
However, the process involves the complication of the cubic eigenvalue problems, and the 
results are usually far from transparent for a clear extraction of the physical implications. 

In this work we study this problem from another angle. It is well-known that, when 
neutrinos propagate through matter, the latter contributes an induced mass to the neutrinos. 
Similar to the case of RGE, we may write down, as a function of the induced mass, a set of 
evolution equations for the neutrino parameters. It turns out that, with the initial conditions 
given by Wo in Eq. (jHJ), we can find simple, approximate, solutions to these equations, as we 
will detail in Sec. IV. These results were summarized in a previous publication 

To derive these equations, we start from the effective Hamiltonian for neutrino propaga- 
tion in matter 

H eff = H/2E, (9) 
where H is given, in the flavor basis, by 



H 
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mi 



Vn 



m 



r r , T 



mi 



(A 



\ 







(10) 



where mi, m 2 , and m 3 are the neutrino 
vacuum, E is the neutrino energy, and the induced mass A 
The matrix H can be diagonalized, 



masses in vacuum, Vq is the mixing matrix in 
V2G F n e E. 



H = VDV ] = V 



(Di 



\ 



Do 



D, 



(11) 



where Di = M 2 is the squared mass in matter. To study how the elements of V evolve in 
matter, one may start with dH/dA, which leads to 



V^iVDVnV 



\V n \ 2 V 12 V& V 13 V^ 
V n V? 2 \V 12 \ 2 VuVi 



13 "12 



(12) 
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1 1 1 n n \ 


l/(i>>2 - ^3 J 


l/(i>>3 - ^l) 


dx\/dA 


xix 2 - 2xiy 2 + 2/12/2 


-xix 2 + xix 3 + 2/12/2 - 2/12/3 


-X1X3 + 2xi2/3 - 2/12/3 


ax2/dA 


-x x x 2 + 2x 2 yi - 2/12/2 


x 2 x 3 - 2x 2 2/ 3 + 2/22/3 


XiX 2 - x 2 x 3 - 2/12/2 + 2/22/3 


dx 3 /dA 


-X1X3 + X 2 X 3 + 2/12/3 - 2/22/3 


-x 2 x 3 + 2x 3 2/ 2 - 2/22/3 


xix 3 - 2x32/i + 2/12/3 


dyi/dA 


-xix 2 + 2x 2 yi - 2/12/2 


-xix 2 + xix 3 + 2/12/2 - 2/12/3 


X1X3 - 2x32/i + 2/12/3 


dy 2 /dA 


x x x 2 - 2x x y 2 + 2/iy 2 


-x 2 x 3 + 2x 3 2/ 2 - 2/22/3 


XiX 2 - x 2 x 3 - 2/12/2 + 2/22/3 


dys/dA 


-xix 3 + x 2 x 3 + 2/12/3 - 2/22/3 


x 2 x 3 - 2x 2 2/ 3 + 2/22/3 


-X1X3 + 2xi2/3 - 2/12/3 


d(ln J)/dA 


-xi + x 2 + 2/i - 2/2 


-x 2 + x 3 + 2/2 - 2/3 


ail - x 3 - yi + 2/3 



TABLE I: dxi/dA, dyi/dA, and d(ln J)/dA are expressed as sums of terms in l/(-Di —D 2 ), 1/(D 2 — 
D 3 ), and 1/(D 3 -Di). 



Taking the diagonal terms of Eq. (Tl2|) , we find 

dDi 



dA 

The off-diagonal terms yield 



Vur = Xi- Vi , (2 = 1,2,3). (13) 



The diagonal elements {(dV^ /dA)V]u are not constrained by Eq. (1121) . Fortunately, it is 
rephasing dependent [l0| , and we can set it to vanish by a proper choice of the phase. This 
means that, when we multiply Eq. ( 114)) by (V')kj, and sum over k ^ i on the right hand 
side, we may sum over all k- values on the left. The result is 

dVij x VikVij 



E Tf^TT^- (is) 



dA ±Z Dj - D k 



Note that the dependence is only on the mass differences, (Dj — D k ), in accordance with 
the invariance of Vij if H — > H + constant. 

While Eq. ( jl5p is valid only with a particular choice of phase, this rephasing ambiguity 
is removed if one uses it to compute rephasing invariant quantities, e.g., 

dVi23 d dxi .dJ 

-lA=dA {VM ^ ) = lA- t dA- (16) 

After some algebra, separating the real and imaginary parts, in addition to using different 
T- k s, we obtain the evolution equations for all (xi,yA and dlnJ/dA, which are collected in 
Table I. 

Note that, since Eq. ( fl5|) can be obtained from Eq.(3.6) (Ref. [lfj]) in appropriate limits, 
the entries in Table I can be identified with those in Table II of Ref. [To| . Indeed, it can be 
verified that d Xi /dA = J2[(Bi) 2r -(Bi) 3r ]/(D s -D t ), dy./dA = Ep^-^KJA^-A), 
where the sum is over cyclicly permuted (r,s,t) = (1,2,3), and B^B'A are obtained from 
Bi(B'j) in Table II of Ref. [Toj by exchanging X2 -H- x 3 , since V -H- under the usual 
conventions in going from quarks to neutrinos. 
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The symmetric form of these equations allows us to find readily the result: 



^ lnfjp! - D 2 ){D 2 - D 3 )(D 3 - £>i)] = 0, (17) 
i.e., the product [J{D\ — D 2 )(D 2 — D 3 ){D 3 — Di)] is a constant as A changes, a well known 



result derived with different methods [H 
From Table I, we find 



~ M* 2 - 2/2) = + (19) 



We see that there is another "matter invariant" : 

d r J 2 



-yi) = 


x 2 - 
D x - 


"2/2 

-D 2 


x 3 - 
D 3 - 


"2/3 

-Dx 




Xi 


-yi 


x 3 
D 2 


-2/3 


Dx 


-D 2 


-D 3 


2/3) = 


x 2 


-2/2 


D 3 


-2/i 


D 2 


-D 3 


-Dx 



dA (xx - yx)(x 2 - y 2 )(x 3 - y 3 ) ' 

Or, 

t2//|t/ |2|t/ |2|t/ 12 



(20) 



0. (21) 



^ /(l^ii f |^12| l^isl ) = constant. (22) 
When we use the "standard parametrization" , it is seen that J 2 /{\Vu | 2 |^i2| 2 |^ // i3| 2 ) = 



121 



S^,S 23 C 23 , i.e., S^sm29 23 is independent of A, a result obtained earlier 1J 

The evolution equations for (x, y) also have a structure akin to that of the fixed point of 
single variable equations. It can be verified that, if + yi — (i — 1, 2, 3), then 

A( Xi + % .) = , j = (1,2,3). (23) 
This result is understandable since the conditions Xi + yi = are equivalent to W 2 i = Wj. 



which, in turn, imply that the effective Hamiltonian H has a \i — r exchange symmetry [13 
This symmetry is clearly independent of A in Eq. ffTU]) . resulting in Eq. ( l2"B"j) . Note also 
that there are actually only two independent constraints in x,- t + = 0. Given any two of 
them, say for i = 1,2, we can use Eq.(4) to derive x 3 + y 3 = 0. Thus, the set of evolution 
equations has a "fixed surface" (in the four-dimensional parameter space), points on the 
surface defined by Xi + y, = stay on it as A varies. 



IV. APPROXIMATE SOLUTIONS 

While analytical solutions to the equations in Table I are not available, as we will see, 
given the known physical parameters, one can exploit certain characteristic properties thereof 
to arrive at simple, but fairly accurate, solutions to these equations. Before we do that, it is 
instructive to first study the two flavor problem, which can be compared to the traditional 
approach, since exact solutions can be obtained in both cases. 
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A. Two-flavor problem 



For two flavors, we have 

dH_flO 
dA ~ \ 

with the familiar diagonalization matrix 



(24) 



/ oos9 sm9\ 

[-sine cos6 J ' K } 

so that x = V11V22 = cos 2 9, y = V12V21 = — sin 2 9, x — y = 1. The evolution equations are 

dD . . 

dA = + 



(26) 



dx 2xy dy 
dA = 1T = dA' 

where D = D 2 — D\. It follows that 

±(xyD>) = 0, (27) 

Jj[D(x + y)} = -(x + y) 2 + Axy = -1. (28) 
Eq. (p7|) is the familiar result 

D 2 sin 2 2^ = D 2 sin 2 20 o . (29) 

Eq. (J2ED gives 

D(x + y) = -A + D (x + y) , (30) 

and thus 

D 2 = [(A - D cos 2# ) 2 + Dl sin 2 29 ], (31) 

which is the well-known resonance formula with D = min. at A = Dqcos29q. 

These results show that the use of evolution equations is equivalent to the traditional 
method, that of finding directly the eigenvalues of the effective Hamiltonian. We now turn 
to the case of three flavors. 



B. Three-flavor problem 

Experimentally, it is known that 5$ = m\ — m\ = 7.6 x 10~ 5 eV 2 , A = m\ — m\ = 
2.4 x 10~ 3 eV 2 , so that 5q/Aq 1 (We assume the "normal" ordering of neutrino masses. 
The "inverted" case can be similarly treated). Note that these values are relevant to long 
baseline experiments since A = V2G F n e E ~ (7.6 x 10- 5 eV 2 )(E/GeV)(p/ gem" 3 ). 

Since 60 A , we expect that the three-flavor problem can be approximated by a pair of 
well separated two-flavor problems [13]. Indeed, the structure of the differential equations in 
Table I shows that the variables (xi, yj) evolve slowly as a function of A except for two regions, 
where Di w D 2 and D 2 ~ -D 3 , corresponding to the two resonance regions. More precisely, 
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let us denote by {A , Ai, A i: A h , A d ) the values of A in vacuum {A = 0), at the lower 
resonance (Ai, [d(Di — D 2 )j dA\^ = 0), the intermediate range {Ai), the higher resonance 
{Ah), [d{D 2 — D 3 )/dA]A h = 0), and for dense medium {Ad). Rapid evolution for only 
occurs for A w A\ and A Ah. Near the lower resonance region, {D 2 — D\) -C (-D3 — D\) or 
(Z?3 — D2). For the higher resonance region, {D3 — D 2 ) <C {D$ — D\) or {D 2 — D{). Thus, for 
these two regions, we need only to keep terms oc 1/ {D 2 — D\) and 1/{D 3 — D 2 ), respectively 
This approximation is generally valid to O(10~ 2 ). We now turn to a detailed analysis. 

For < A < Ai, in the neighborhood of Ai, we keep terms oc \/{D\ — D 2 ) in Table I. 
Let's concentrate on the variables Xi — yi — \Vu\ 2 and use Eq. ( TT5]) and Eqs. (IT5|) - (j2D|) . We 
define 



Then, 





X 


= Xi - 


2/1, 




Y 


= x 2 - 


V2, 




Z 


= x 3 - 


y-3, 




5 


= D 2 - 


-D x . 


dX 




2XY 




~d~A 




5 




dZ 
dA 




0. 




d5 
dA 






Y). 



(32) 



dY 
dA'' 



(33) 

These equations are identical to those for the two flavor problem, Eq. (12 6p . with X — > x, 
Y — > —y, and 5 — >■ D. Also, in place of x — y = 1, we have 

X + Y = Pl , (34) 

where pi is a constant (since d{X + Y)jdA = 0), and pi = 1 — Z, Z = \Vi 3 \ 2 = constant. 
The solutions are 

XY5 2 = X Y 5 2 , 
{X-Y)5 = -p 2 A + gi 6 , 

qi = X - Y . (35) 

Explicitly, we have 

5 2 = pfA 2 -2 qi 6 A + 5 2 , 
X = l-[ Pl -(p 2 A- qi 5 )/6}, 

Y = l -[p l + (p 2 A- qi 5 Q )/5}. (36) 

Thus, {5, X, Y) exhibit the classic resonance behavior, with the resonance location at the 
minimum of 5: 

A t = {%)S . (37) 
Pi 



S 



Substituting in the vacuum input values, X = 2/3, Y = 1/3, A\ = S /3. The width of the 
resonance is 

(SA) l = (l-qf/p^ 2 (5 / Pl ). (38) 

For the physical PMNS matrix, (5A)i ~ 5q. This means that the intermediate A value, Ai, 
already starts at A > (2 - 3)6 . For A ~ A i: 5 ->■ A, X ->■ 0, K -> 1, with ^ = 1. 
Turning to the higher resonance, we define 



The evolution equations are 



with the solutions 



dX 

~d~I 

dY 

dA 
dA 

dA 



A = D 3 - D 2 . 



0. 



2FZ 



(39) 



(40) 



YZA 2 = Y Z A 2 , 
(Y-Z)A = - P 2 h A + q h A , 



(41) 



or, 



A 2 
Y 



PhA 2 - 2q h A A + Ag, 
^[ Ph -(p 2 h A-q h A )/A], 

^[ Ph +(p 2 A~q h A )/A}, 



(42) 



where ph = Y + Z , q^ = Y — Z . Here, the values Fo and Z are taken at A = Ai ^> So, 
so that Yq = 1, Zq = \Vis\ 2 = 0, from the solutions for the lower resonance. The position of 
the higher resonance is at 

A h = (^)A = A . (43) 



Its width is 



PI 



(SA) h = (1 - ql/pD^iAo/ph) = 2Z A « A . 



(44) 



The above analyses show that the two-flavor approximation yields simple solutions to 
the mixing parameters |Vij| 2 , for all A values. However, the vacuum mixing, given by Wo 
in Eq. dHJ), has another important feature, namely, (Wo)2« — {Wo)^, or Xio + yio — 0. This 
feature, according to Eq. ([23"]) . is preserved by the evolution equations and so W 2 i — W 3i , or 
+ y i = 0, for all A. Thus, with the known solutions for Wu from above, all elements of 



li 

W are determined by unitarity. Explicit solutions for W or (xi,yi) were presented in Ref. 
obtained by using both approximations. 

We may divide the full range of A values into a low- A and a high-v4 regions. The former 
covers the range from A = to a value below the higher resonance region, while the later 
starts from beyond the lower resonance region and ends at A = oo. In these regions, the 
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FIG. 1: Numerical (solid) and approximate (dot-dashed) solutions for (a) all D 3 {A), D2(A), and 
Di(A), and (b) the enlarged plot of D 2 {A) and D X (A) in < A/5 < 10. We have used 5 /A = 
1/32. 

Di/6 Di/6 




evolution equations are dominated by contributions from pole terms, 1/(D\ — D 2 ) for low- 
A and \j(D 2 — D 3 ) for high-A. The exact demarcation between low- A and high-A is not 
important, since in the intermediate region contributions from either pole are small, and 
there can be considerable overlap between low- A and high-A, corresponding to the large 
range of Aj. 

It should be noted that "pole dominance", which was used to go from Eqs. f fT8ll20|) to 
Eqs. (155j) and (|40p . is an excellent approximation in this situation. This is because the 
terms dropped are doubly suppressed, first by the large denominators, and then by the 
small numerators ((x 3 — y 3 ) <C 1 throughout the low- A region and (x\ — y\) <C 1 for high- 
A). We should also emphasize that Eqs. (133|) and ( 140]) are derived independently of the 
approximations Xi + yi ~ 0. While in our earlier work [9|, they were used to arrive at similar 
equations of evolution. 

The "pole dominance" approximation, less accurately, may also be used for other vari- 
ables, such as J 2 . This gives rise to relations, which we shall dub as "partial matter invari- 
ants" , valid only for either the low- A or the high-A regions. Thus, from Eq. (155]) and using 
Table I, 



(Di - D 2 y\V 11 \ 2 \V 12 \ 2 = constant, (low-A) 

J 2 (D 1 — D 2 ) 2 = constant. (45) 



Similarly, 



(D 2 - D 3 ) 2 \V 12 \ 2 \V 13 \ 2 = constant, (high-A) 

J 2 (D 2 -D 3 ) 2 = constant. (46) 

These "partial matter invariants" are useful in understanding some detail properties of the 
parameters. E.g., for A ~ A4, \Vn\ 2 = (2/9)(5q/A) 2 . The behavior of J 2 is also clarified, as 
we will see in the discussion on Fig. 4. 

In summary, the solution to the three flavor problem can be made simple by dividing the 
full range of A into a low-A and a high-A regions, In the low-A region, the evolution equations 
are dominated by pole terms oc i/(D% — D 2 ), and the solution centers around the lower 
resonance. For the high-A region, correspondingly, pole terms oc 1/{D 2 — D 3 ) dominate, 
and the solution can be characterized by the higher resonance. The mixing parameters 
change appreciably only in two regions: 1) lower resonance, [A/ — (SA)i\ < A < [A/ + (<L4)j], 
At = (Jo/3, (5A) t = 5 ; 2) higher resonance, [A h - {5A) h } < A < [A h + [5A) h ], A h = A , 
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(5A) h S 2|y 13 | 2 A < A . The solutions for X(|y n | 2 ), K(|y 12 | 2 ), and Z(\V 13 \ 2 ) are: 1) 



A : 

X 
X 
X 



0, X 2/3, F = 1/3, Z 
F = 1/2, Z = e 2 = 0; 3) A 



0, y 1, Z = e 2 ; 4) A 
Y = 0, Z = 1. 



' e 2 < 1, which are the given vacuum values; 2) A = A;, 
Aj, Aj covers roughly the range, 2<5 < Aj < A (l — 2e 2 ), 
X = 0, y Z = 1/2; 5) A = A d , with A d > A (l + 2e 2 ), 



When we incorporate the other feature of the vacuum PMNS matrix, that (Wo)2i — 
(W / o)3i, which is preserved by the evolution equations, the result is that Wn — W 3 i, for all 
A. Given Wu from above, the matrix W is then completely determined by unitarity. 

Our results can be put together by giving the matrices W at A = (A , A[,Ai,A h , A d ): 



Wn 



w. 



f 2/3 1/3 \ 

1/6 1/3 1/2 
^1/6 1/3 1/2 ) 

' 1 \ 
1/2 1/2 , 

V 1/2 1/2 ) 

' l\ 
1/2 1/2 . 

V 1/2 1/2 Oj 



' 1/2 1/2 \ 
i 1/4 1/4 1/2 
V 1/4 1/4 1/2 ) 

' 1/2 1/2 \ 

1/2 1/4 1/4 , 
^ 1/2 1/4 1/4 ) 



(47) 



As a group, these matrices exhibit the remarkable simplicity of the PMNS matrix as A 
varies from to oo. Note that all of the matrices have at least one zero, Wn = 0, implying 
Xi = yj = 0. 

The other feature, as mentioned before, is that they have equal elements in their second 
and third rows, Wn = W^. This means that the W matrix is completely fixed by its first 
row, Wu- These elements, in turn, control dDi/dA, Eq. ([IB"]) . Thus, the progression of W as 
a function of A can be read off from the plot of Di(A), which is given in Fig. 1. Generally, 
the accuracy of the entries in Eq. ( 14T|) is of the order 10~ 2 . Note also that A« covers a 
rather large range, 28q ^ Aj ^ A (l — 2e 2 ). The validity of Eq. (T4T|) will be confirmed by 
numerical integrations, to be given in the next section. Where applicable, they also agree 
with numerical results in the literature jsj, after the proper change of variables is carried 
out. 



V. NUMERICAL SOLUTIONS 

It is straightforward to numerically integrate the evolution equations for (x,y). To do 
this we first obtain the vacuum expressions for the (x,y) parameters from Eqs. (6) and (8): 

x w = ~(2 - 3/3 - 2e 2 ), y w = ~(-2 - 3/3 + 2e 2 ), 
o o 

x 20 = J(l - 3/3 - e 2 ), y 20 = -3/3 + e 2 ), 

o o 

x 30 = ^((3 + e 2 ), 2/30 = I(/3- e 2 ), (48) 

where £ = rj = is chosen and the terms in 0((3e 2 ) are ignored. In addition, we choose the 
initial values e = 0.17, (3 = 0.02, corresponding to the experimental bounds | V^ 3 1 2 < 0.03 [BJ 
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FIG. 2: The numerical (solid) and approximate (dot-dashed) solutions for x±(A), X2(A), and xs(A). 
Note that yi(A) ~ -Xi{A). 
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and an assumed CP violation phase cos0 =1/4. With the input hierarchy <5o/Ao = 1/32, 
the numerical results for Di and that for (xi,yi) are then compared with the approximate 
solutions obtained earlier (Eqs. f )36|) . f )42|) ) in Fig. 1 and Fig. 2, respectively. The agreements 
are quite good. 

Given the possible normal or inverted mass hierarchies, Fig. 3 summaries the evolution 
of the mixing parameters for both the v and the v sectors. Note that the parameters (x~i, f/i) 
for the v sector can be obtained by replacing A with —A and V with V* in that for the 
v sector. With the establishment of these "basic solutions" for the mixing parameters, we 
may easily study other physical quantities that are relevant to practical calculations for the 
neutrino propagation in matter. For illustration purpose, we plot some of the quantities 
numerically under the normal hierarchy in the following. The corresponding solutions under 
the inverted hierarchy can be manipulated likewise. 

The evolution of J 2 = Xix 2 x 3 — 2/12/22/3 {J 2 — Xix 2 x 3 — 2/12/22/3) m matter is shown in Fig. 
4. Compared to its vacuum value, it is seen that, except for some enhancement for J 2 near 
A = Ai and A = Ah, the general trend is for it to decrease with A. We note that near Ai, in 
the l/(Dx — D 2 ) dominance approximation, J 2 [D\ — D 2 ) 2 ~ constant, so while (D 2 — Di) 
goes through a dip, J 2 has a bump near Ai, after which J 2 / Jq w (<5 / A) 2 , for <5 A < A . 
A similar behaviour occurs near Ah, with J 2 (D 2 — D3) 2 w constant, although the effects 
are hardly noticeable. Similar numerical results were also reached by solving directly the 
eigenvalue problem [8j. In addition, Fig. 5 shows the evolution of the first two rows of Wij 
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FIG. 3: The mixing parameters in matter under the normal hierarchy (left column) and the inverted 
hierarchy (right column), for both the v sector (solid, Xj) and the v sector (dashed, Xj). Note that 
Xi + yi ~ (xi + yi ~ 0) for all A, and X, = x, in vacuum. 
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FIG. 4: The evolution of J 2 = X1X2X3 — 2/12/22/3 and J 2 = X1X2X3 — 2/12/22/3, obtained from the 
numerical solutions for v (solid) and v (dashed) sectors. 
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in matter. We do not present the plots of Wz% since they are almost indistinguishable from 
those of W2i- The patterns of in Eq. ( H7|) are clearly seen from the plots. Furthermore, 
the mixing angles of the standard parametrization: sin 2 12 , sin 2 #23, and sin 2 #13, are related 
to the (x, y) parameters: 

sin 2 fli2 = l/(l + ^^), (49) 

X 2 ~ 2/2 
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FIG. 5: Wij for both the v (solid lines) and v (dashed lines) sectors are plotted as functions of 
A/5q. Note that Wn = W$i, and only the first two rows of Wij are shown in the figure. The 
patterns of Wij can be readily deduced from Fig. 3. 



Wu W 2 i 




FIG. 6: The mixing angles for the v (left) and the v (right) sectors in the standard parametrization: 
sin 2 #12 (solid), sin 2 #23 (dot-dashed), and sin 2 #13 (dashed), are plotted using the (x,y) values in 
Fig. 3. Note that #23 for the v sector remains constant (#23 = vr/4) for A < A^ (higher resonance) 
and deviates slightly from tt/4 for A > A^. 




X2 - yi 



sin 2 0i3 = x 3 - y 3 . (51) 

The numerical results for both the v and v sectors are shown in Fig. 6. These plots are 
in agreement with those in the literature [§J. For 6 12 and # 13 , note the characteristic step- 
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function resonance behaviors near Ai and Ah. Also, #23 — 7r /4 is a reflection of W23 — W33, 
for all A. The phase angle (0) is not plotted since it also remains constant due to the 
invar iance of S^sin2# 2 3- 

VI. CONCLUSION 

Understanding the propagation of neutrinos through matter is one of the core problems 
in neutrino physics. In a medium of constant density, it is well-known that the electron 
neutrino acquires an induced mass which alters both the eigenvalues and the mixing matrix 
of the neutrinos. Traditionally, one studies directly the eigenvalue problem of the effective 
Hamiltonian. The neutrino parameters are expressed as complicated formulas in terms of the 
induced mass and their values in vacuum. One then resorts to numerical plots by assuming 
specific values for these partially known parameters. The drawback of this method is the 
lack of insights into the nature of the solutions, and it is not easy to gain an overview of the 
mixing as a function of the induced mass. 

In this paper we try a different approach, by finding the evolution equations of the 
neutrino parameters as a function of the induced mass. The resulting equations, when 
written in terms of a rephasing invariant parametrization, turn out to be manageable and 
we are able to find simple, approximate, solutions with the help of two important features of 
the vacuum neutrino parameters. 1) The two measured mass differences are widely separated 
so that the two-flavor resonance approximation becomes applicable. 2) The vacuum PMNS 
matrix has an approximate \l — t symmetry, which is preserved by the set of evolution 
equations. The result is summarized in Eq. (I47p . showing the striking simplicity of the 
neutrino mixing matrix as a function of A. Approximate solutions for the parameters are 
explicitly given in Eqs. (136 j) and (142 p . The evolution equations also facilitate the derivation 
of "matter invariants", given in Eq. ( ITT)) and (122)) . In addition, there are also "partial 
matter invariants", Eqs. (T4"5)) and (|46j) . These are useful in obtaining properties of the 
various parameters without performing detailed calculations. 

Based on the incomplete measurements that exist for the vacuum parameters, our analy- 
ses show that those in matter, to a good approximation, can already be determined. We hope 
that these results will be helpful in the exploration of the physics of neutrino propagation 
in matter. 
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Appendix A: Neutrino transition probabilities in (x, y) parameters 

As the neutrinos travel through a baseline L in matter of constant density, the flavor 
transition probability is given by 

P{y a Uf} ) = S a0 - A^ReiV^V:^) sin 2 (Ai) 

j>i 
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+ 2 £ sin(2A,), (Al) 

where = (Dj — Dj)L/AE. For a ^ we obtain the explicit expression, 

PK->^)= - 4[Re(V al V^V: 2 V^ 2 )sm\D u ) + ReiV^V^V^sm^D^) 

+ Re(V a2 V; 2 V: 3 Vp 3 )sm 2 (D 23 )] 

+ 2[Im(V al V; i V: 2 Vp 2 ) sin(2L> 12 ) + Im(V al V p \V: 3 V^ 3 ) sin(2Z> 13 ) 

+ Im(V a2 V; 2 V: 3 Vp 3 )sm(2D 23 )]. (A2) 

For a specific process, e.g., P{y^ —> v e ), we have 

Re(V 2 iV 12 V 22 V^) = x 2 x 3 + xiy 2 - y x y 2 - y 2 y 3 = Fg, 
Re{V2iV 13 V 2 \V? x ) = -x x x 3 - x 2 x 3 + x 3 y x + y 2 y 3 = Fj£ , 
Re{V 22 Vi 3 V 23 V* 2 ) = xix 3 + x 2 y 3 - 2/12/3 - y 2 y 3 = Fg 2 e , (A3) 

and the probability, 

P{y„ 4[F 2 f sin 2 (£> 21 ) + i^ 6 sin 2 (L> 31 ) + sin 2 (L> 32 )] 

+ 8Jsin(D 21 )sin(Z) 31 )sin(D 32 ), (A4) 

where Jw[V^ V^j V^- VJJ = J X) 7 ,fe e ap~/ e ijk h as been used in reducing the sum of the imaginary- 
parts in Eq. (1A2j) . In addition, the probability for the T-conjugate process takes the form 



P{ve -+vy)= - 4[F 2 T sin 2 (i) 21 ) + F 3 T sin 2 (L> 31 ) + F 3 T sin 2 (L> 32 )] 

- 8Jsm(D 21 )sin(D 31 )sm(D 32 ), (A5) 



where 



F 2i = ~x\x 2 - xix 3 + x x y 2 + 2/12/3, 

F 3i = x i x 2 + x 3 yi - ym - ym, 

F 32 = -xix 2 - x 2 x 3 + x 2 y 3 + yiy 2 . (A6) 

We may verify the relation Fj^ e = F^ using Eq.(4). 

The explicit probabilities for other processes can be derived following the same procedure: 

P(v e -)■ v T ) = - 4[(xix 3 + x 2 y x - 2/12/2 - J/12/3) sin 2 (A>i) 
+ - x x x 3 + x x y 3 + 2/12/2) sin 2 (£) 3 i) 

+ (xix 2 + ^32/2 - 2/12/2 - 2/22/3) sin 2 (£> 32 )] 

+ 8Jsin(i5 21 )sin(D 31 )sin(D 32 ) (A7) 



P{v T ->v e )= - 4[(-xix 2 - x 2 x 3 + x 2 yi + 2/22/3) sin 2 (D 2 i) 

+ (x 2 x 3 + Xi2/ 3 -2/12/3 -2/22/3) sin 2 (D31) 
+ - x 2 x 3 + x 3 y 2 + 2/12/3) sin 2 (£> 32 )] 

- 8Jsinp 21 )sinp3i)siiip32) (A8) 
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P(fp v T ) = - 4[(-xiX 3 - x 2 x 3 + x 3 y 3 + y x y 2 ) sin 2 (-D 2 i) 
+ (x 1 x 3 + x 2 y2 - 2/12/2 -2/22/3) sin 2 (D31) 
+ (-xix 2 - XiX 3 + xij/i + 2/22/3) sm 2 (D 32 )] 
- 8Jsm(D 21 )sm(D 31 )sm(D 32 ) (A9) 



P(i/ T -^v li )= - 4[(xiz 2 + ^32/3 - 2/12/3 - 2/22/3) sin 2 (D 2 i) 

+ (-X1X2 - x 2 x 3 + x 2 y 2 + 2/12/3) sin 2 (P 3 i) 

+ (0:2X3 + xxyi - 2/12/2 - 2/12/3) sin 2 (D 32 )] 

+ 8Jsin(i52i)sin(i53i)sin(D32) (A10) 

We may also write down the expressions for the v sector, P{p a — > up), by replacing the 
parameters for the v sector with that for the v sector: x — >■ x, y — )■ 27, -Dy — >■ -Dij, and thus 
— >■ Pj" , J — J. As an example, the probability Piv^ — >■ z? e ) is given by 

Pfo * e ) = - 4[F 2 T sin 2 (,D 21 ) + F3T sin 2 (An) + P32 sin 2 (D 32 )] 

- 8Jsm( J D2i)sm( J D 3 i)sm(L>32), (All) 

where 

P2T = ^2X3 + x x y 2 - 2/12/2 - 2/22/3, 
^3i £ = -X1X3 - x 2 x 3 + X32/1 + 2/22/3, 

P32 = X1X3 + x 2 y 3 - 2/12/3 - 2/22/3- (A12) 

The evolution equations and the analytic, approximate, solutions for (xj, j7j) can be obtained 
following the same method outlined in this work. The numerical solutions for (xj,f/j) are 
shown in Fig. 3. Note that 1) = F^ = F^ a = F^ a in vacuum; and 2) the functions 

F°f ', Fij 01 , F°j v and F^ a can take varied forms in terms of (x, y) and (x,y) since there are 
different ways |2| of reducing Re(Vc^VgjV*jVpj). 
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